

//u1 = sqrt((ZA - Z1 + PA / (1000 * 9.81) + uA * uA) / (1 + (0.02 * le1) / d));

